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We investigate the initial state of the inflationary universe. In our recent publications, we showed that request- 
ing the gauge invariance in the local observable universe to the initial state guarantees the infrared (IR) regularity 
of loop corrections in a general single clock inflation. Following this study, in this paper, we show that choosing 
the Euclidean vacuum ensures the gauge invariance in the local universe and hence the IR regularity of loop cor- 
rections. It has been suggested that loop corrections to inflationary perturbations may yield the secular growth, 
which can lead to the break down of the perturbative analysis in an extremely long term inflation. The absence 
of the secular growth has been claimed by picking up only the IR contributions, which we think is incomplete 
because the non-IR modes which are comparable to or smaller than the Hubble scale potentially can contribute 
to the secular growth. We prove the absence of the secular growth without neglecting these non-IR modes to a 
certain order in the perturbative expansion. We also discuss how the regularity of the n-point functions for the 
genuinely gauge invariant variable constrains the initial states of the inflationary universe. These results apply 
in a fully general single field model of inflation. 
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I. INTRODUCTION 
A. Motivation and the current status of IR issues 



Initial states of the observable universe. How our universe began? This is one of the biggest question in cosmology. 
The observation of the cosmic microwave background tells us about the cosmological perturbation at the last scattering, which 
realistic scenarios of the early universe should explain. If inflation took place preceding the big-bang nucleosynthesis, the 
OO " quantum fluctuation of the inflaton can generate the seed of the cosmological fluctuation which is consistent with the scale- 
QO invariant spectrum at large scales. Therefore, in the context of the inflationary scenario, which is currently the most successful 
scenario of the early universe, the period of inflation is the earliest part of our observable universe. In this series of papers, 
we pursue the question; "What can we claim about the initial quantum state of the observable universe if we require that the 
theoretical prediction should be stable against the infrared (IR) loop contribution?" The adiabatic vacuum is widely accepted to 
be the most natural vacuum at least for a free field theory, since it mimics the vacuum of the flat spacetime in the ultraviolet (UV) 
limit. However, in a number of publications lflj-l23ll. it has been suggested that the adiabatic vacuum may not be stable against 
. . ' the IR contributions in the presence of non-linear interactions. 

Non-locality of the action and IR divergence problem. When we assume that the free field has the scale invariant spectrum 
in the IR limit, a naive consideration can easily lead to the IR divergence due to loop corrections. Here we illustrate how the IR 
divergence can appear from the loop corrections of the curvature perturbation in single field models. Choosing the time slicing 
on which the inflaton field is homogeneous, we can express the action in terms of the unique dynamical degrees of freedom (, 
the curvature perturbation, and the Lagrange multipliers TV and TVi, the lapse function and the shift vector. The Hamiltonian 
and the momentum constraint equations relate the dynamical variable £ to the multipliers TV and TV,. As is explicitly shown in 
various papers, for instance in Ref. ||24| - |2o1| . these constraint equations are elliptic-type equations, and schematically written as 

d 2 N = f[(], d 2 N i = f i [C\, (1.1) 

where d 2 denotes the spatial Laplacian. By requesting the regularity at the spatial infinity, the boundary conditions of these 
elliptic-type equations are uniquely fixed. Substituting the expressions of TV and TV into the action, we obtain 

S = J d 4 x£[C, TV, TV,] = J d 4 x£[C, d~ 2 f[C], d- 2 fi[Q] , (1.2) 



'Electronic address: tanaka'at'yukawa.kyoto-u. ac.jp 
tElectronic address: yurakawa'at'ffn. ub.es 



2 



and hence the evolution of ( is described by the above non-local action. Here the inverse Laplacian <9~ 2 is usually supposed to 
be defined as multiplying the inverse of the eigenvalue of the Laplacian operator by using the harmonic decomposition. When 
we evaluate the loop corrections to the n-point functions expanding them in terms of the interaction picture field (/, we need to 
evaluate the expectation values such as 

<C?>, <C/d- 2 C/>, •••• (1.3) 

Inserting the scale invariant spectrum into ((f) leads to the logarithmic divergence as ((f) oc J d 3 k/k 3 . The second expression 
of Eq. ( 11.31 ). which may arise as a consequence of the operation of d~ 2 , is more singular as ((/9~ 2 (j) oc J d 3 k/k 5 , which 
diverges quadratically. The presence of non-local interactions enhances the long range correlations, and hence the singular 
behaviour in the IR. When we introduce the IR cutoff, say at the Hubble scale at a particular time to, the variance ((f) shows 

the logarithmic secular growth as ((f) oc J^ H H dk/k ~ loga/ao where ao and Hq, respectively, denote the scale factor and 
the Hubble scale at t = to- If the IR divergence exists, the loop corrections, which are suppressed by an extra power of the 
amplitude of the power spectrum (H/M p \) 2 , may dominate in case inflation continues sufficiently long, leading to the break 
down of perturbation. 

The dilatation symmetry as a necessary ingredient for IR regularity. The regularization of the IR contributions has been 
discussed in a number of publications I25K38I1 . The important aspect in discussing the long wavelength mode of ( is the dilatation 
symmetry of the system. As is expected from the fact that the spatial metric is given in the form a 2 e 2< *dx 2 , a constant shift of 
the dynamical variable ( can be absorbed by the overall rescaling of the spatial coordinates. Hence, the action for ( preserves 
the dilatation symmetry: 

x l -> e~ s x i , ((t, x ) -> CO, e" s £c) - s , (1.4) 

where s is a constant parameter. (There are a number of literatures where this dilatation symmetry is addressed. See for 
instance, Refs. IT391 l40ll and the references therein.) One may naively expect that we can absorb the IR divergent contribution 
of ( using this constant shift. As an example, we set the parameter s to the spatial average of the curvature perturbation within 
the Hubble patch at to, ((to), where the size of the Hubble patch in comoving coordinates is given by 1/(oqHq). Then, the 

logarithmically divergent two-point function ((f) seems to be replaced with (((j — (~/) 2 ) oc f£ H dk/k, which is finite but 
still grows logarithmically in time. One may think that if the system is described in such a way that the symmetry under the 
time dependent dilatation transformation is manifest, setting s(t) to the time dependent spatial average in the Hubble patch, the 
logarithmic growth of ((f) might be eliminated. However, the reduced action written in terms of ( dl.2t does not preserve the 
invariance under the dilatation transformation with the time dependent parameter s(t). For example, in the recent literature ll40ll . 
the authors showed that when we consider the whole universe with the infinite spatial volume, the dilatation transformation 
should be time independent to preserve the action invariant. In addition, the two-point function with d~ 2 cannot be regularized 
by considering the dilatation symmetry alone. This quick consideration tells us that the presence of the dilatation symmetry 
of the system may play an important role in the regularization of the IR contributions but is not sufficient to guarantee the IR 
regularity and the absence of the secular growth. 

Residual gauge degrees of freedom in the local universe. A missing piece in the above discussion is to pay careful attention 
to what are the quantities we can actually observe. Since our observable region is a limited portion of the whole universe, the 
observable fluctuations must be composed of local quantities. Furthermore, as the information that we can access is limited to 
our observable region, there is no reason to request the regularity at the spatial infinity in solving the elliptic constraint equations 
( 11. It . Then, there arise degrees of freedom in choosing the boundary conditions of Eqs. (Il.lt . The degrees of freedom in 
solutions of N and Ni can be understood as the degrees of freedom in choosing coordinates. As we showed in Refs. I25ll26ll . 
these residual coordinate transformations are expressed in terms of homogeneous solutions to the Laplace equation as 

X i -> X i - S (ty - Si h-3m (t) xh ■ ■ ■ ^ + ■ ■ ■ > (1-5) 

m— 1 

where s l j 1 ...j m (t) are symmetric traceless tensors, which satisfy 5^ s l j 1 ---j—j'---j m {t) = 0. Here, we abbreviated the non- 
linear terms in the coordinate transformation. Note that this coordinate transformations include the dilatation transformation 
with the time dependent function s(f). Since the transformations in Eq. (11. 5t are nothing but coordinate transformations, the 
diffeomorphic invariant action S = J d 4 xL[(, N, Ni] should preserve the symmetry under these transformations. Thus, when 
we consider only the local observable region, which is a potion of the whole universe, we find an infinite number of coordinate 
transformations which keep the action invariant. Considering the dilatation transformation in the whole universe is subtle in 
the sense that the transformation diverges at the spatial infinity, even if the parameter s is very small. By contrast, restricted 
to the local region, the magnitude of the coordinate transformations in Eq. ( 11.5b is kept perturbatively small. In this paper, we 
refer to the local observable (spacetime) region as 0. The size of the observable region on each time slicing is supposed to 
be of order \ / a(t)H(t) at least in the far past since the past light cone asymptotes to that size. We should note that, once we 
insert the expressions of N and N{ into the action to obtain the action for the curvature perturbation (, the symmetry under the 



3 



residual coordinates transformation is lost, because specific boundary conditions are chosen for TV and Ni in fixing coordinates. 
To emphasize the distinction between the coordinate transformations associated with the change of the boundary conditions and 
the usual gauge transformation, which keeps the action invariant, we denote the former by the gauge transformation in the italic 
font. 

Removing the residual gauge degrees of freedom. One way to realize the invariance under the gauge transformation is fixing 
the gauge conditions completely. The residual gauge degrees of freedom introduced above can be also removed by employing 
additional gauge conditions, i.e., by fixing the boundary conditions of N and Ni at the boundary of the local region 0. Then, 
we naturally expect that the IR regularity may be explicitly shown by performing the quantization in this local region, since the 
wavelengths that fit within this local region will be bounded by the size of the region. Although the quantization in the local 
region is an interesting approach, it is not so clear how to perform the quantization after removing the residual gauge degrees 
of freedom. One of the difficulties is that even the translation symmetry of the quantum state cannot be easily guaranteed in the 
local system, since it is manifestly broken by introducing the boundary condition at a finite distance. (See also the discussion in 
Ref. S). 

As an alternative way, in Ref. (3, we first set the initial state considering the whole universe, and then we performed the 
residual gauge transformation (II. 5t to fix the coordinates so that the IR contributions are absorbed. Through the transformation 
with s(t) = C(t)> m e curvature perturbation is transmitted as 

C(i, x) -> £(*, e-<®x) - t(t) - C(t, x) - CW + 0(C 2 ) ■ (1-6) 

Here, £(x) is the original curvature perturbation defined in the whole universe and its spatial average over the whole universe is 
set to as in the conventional cosmological perturbation theory. By contrast, £(t, e~^x) — C,{t) is the curvature perturbation 
relevant to the local universe, and its spatial average over the local region S f n is set to 0, where St is a time constant 
surface. In Ref. ll28ll we considered the fluctuation of the inflaton, using the flat gauge, but the same discussion follows also 
for the curvature perturbation £. In the recent publication by Senatore and Zaldarriaga ll38l , the same degrees of freedom in 
choosing coordinates are used in a slightly different way to absorb the IR divergent contributions. If the non-linear terms in the 
residual gauge transformation at the initial time (11.6) did not yield IR divergent contributions, the discussion in Ref. [|28il would 
have proved the absence of IR divergence in general. What was shown there is that once the field operator after the residual 
gauge transformation is guaranteed to be regular at the initial time, its succeeding evolution does not produce IR divergence. 
The heart of the proof is that Ci( x ) is replaced with Ci( x ) — Cf(*) m the expansion of the composite operators in terms of the 
interaction picture field, after the residual gauge transformation, and hence the IR contributions from f j (x) are always canceled 
by those from (i(t). However, the non-linear part of the transformation at the initial time contains whose IR contributions 
logarithmically diverge. The lesson is that it is not straightforward to reformulate the way of quantization so that the IR divergent 
contributions therein are all absorbed by the residual gauge transformation. (The absorption of the IR modes of the curvature 
perturbation was intended in other frameworks such as 5N formalism l3ll l32ll and the semi-classical approach ll33l . ) 

The secular growth. The appearance of IR divergence due to the residual gauge transformation mentioned above might be 
evaded by sending the initial time to the past infinity. This is because the size of the local region E t n in comoving coordinates 
becomes infinitely large in this limit, making the discrepancy between the average in the local region and that in the global 
universe smaller and smaller. Then it might be effectively unnecessary to perform the residual gauge transformation at the initial 
time, although this statement is not very rigorous. We should note that when we send the initial time to the past infinity, it is too 
naive to neglect the non-IR modes which are comparable to or shorter than the Hubble length scale, because all the modes were 
much shorter than the Hubble length scale in the distant past. This makes the issue regarding the secular growth much more 
complicated. For instance, once we include the contributions from the non-IR modes, we cannot use the conservation of in 
the limit k/aH <C 1, where k is the comoving wavenumber of the external leg, relying on the long wavelength approximation 
such as SN formalism. Here, in a simple example, we show that vertex integrations can yield the apparent secular growth 
through the non-linear contributions from the modes at around the Hubble scale. Even if the vertex is confined in the region 0, 
the integration region of each vertex is still infinite in the time direction as J dtd 3 xa 3 (- ■ ■ ) ~ J d(ln a)/H A {- ■ ■ ), which may 
cause the secular growth. Roughly speaking, the integrand (• • • ) will be written in terms of the dimensionless time dependent 
slow roll parameters and the wavenumber of the fields in this vertex k m /aH normalized by the Hubble scale. If we focus on the 
non-linear interaction composed of the modes with k m /aH of order unity, the integrand (• • • ) are expressed only in terms of the 
parameters which are supposed to change very slowly in time and then the contribution from the interaction vertex seems to yield 
the logarithmic growth. This is another origin of the secular growth, which should be distinguished from the one inherited from 
the IR behavior of ((Cr) 2 )- Of course the above argument is too native, but it shows that the absence of the secular growth from 
the vertex integration is rather subtle, requiring more careful treatment about the modes around the Hubble scale. Because of this 
subtlety, introducing the UV cutoff at the length scale longer or equal to the Hubble length scale by hand makes the discussion 
incomplete. In fact, if it were allowed to simply neglect the short wavelength modes, the discussion in Ref. [28] with the initial 
time ti sent to — oo would have given a rough proof of the absence of IR divergence without any limitation to the quantum state 
by sending the initial time to the past infinity, which contradicts our current claim that the quantum state is restricted in order 
to avoid IR divergence. Recently, the absence of the secular growth was claimed relying on the conservation of the curvature 
perturbation in Refs. i37ll38ll . but the aspects mentioned above were not discussed. In addition, even if the conservation of 
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in the limit k/aH <C 1 is proved, the logarithmic enhancement in the form (k/aH) 2 \n(k j aiHi) may give rise, where ai and 
Hi are the scale factor and the Hubble parameter at the initial time. The factor \n(k/a,iHi) can become large to overcome the 
suppression by (k/aH) when we send the initial time to the past infinity. 

B. Summary of upcoming results 

Short summary of the results. In this subsection, we summarize what we will show in this paper. Taking account of the 
current status of IR issues mentioned above, we will establish the following three statements in this paper: 

1. There is an alternative equivalent Hamiltonian that describes the quantum dynamics of our interest and whose interaction 
part is solely composed of the IR irrelevant operators(, which mean the field operators associated with the operations that 
manifestly suppress the IR contribution such as di/aH and d t /H). 

2. The Euclidean vacuum state, which is specified by the regularity when the time coordinates in the n-point functions are 
analytically continued to the imaginary in the complex plane, is physically the same both in the alternative description 
mentioned in item 1, and in the original description. 

3. The n-point functions in the Euclidean vacuum state respect the spatial translation invariance and are regular in the IR. The 
secular growth is absent, even if we include the vertices with non-IR modes, as long as very high order of loop corrections 
are not concerned. 

Below we add a little more detailed explanations about the above three items. 

Gauge issue. In this paper, the quantization and fixing the initial quantum state as a starting point of our discussion is 
performed in the original system which describes the whole universe, where the residual gauge degrees of freedom are left 
unfixed. Then, following Refs. l25U27i[30ll4lll . we introduce a field operator which preserves the invariance under any spatial 
coordinates transformations, including residual gauge transformations. We refer to such an operator as a genuine gauge invariant 
operator. As a representative, we consider a genuine gauge invariant curvature perturbation, 9 R. As long as the expectation values 
of such genuine gauge invariant operators are concerned, we can perform the residual gauge transformation without affecting the 
results of computations. We will show that, using this residual gauge transformation, the boundary conditions of the non-local 
operator d~ 2 in the action can be modified to be regular in the IR. 

Requirement of the gauge invariance in quantum state. To calculate the n-point functions which preserve the invari- 
ance under the residual gauge transformations, the initial state should be also specified in a genuinely gauge invariant manner. 
However, when we perform the quantization considering the whole universe, preserving the residual gauge invariance becomes 
obscure, because these residual gauge degrees of freedom are not present as long as we deal with the whole universe. In our 
previous paper 1127 1, we discovered a correspondence between the IR regularity and the invariance under the residual gauge 
transformations, which will provide an important clue to the guiding principle in choosing the genuinely gauge invariant initial 
state. To discuss this point, aside from the original canonical variables and its conjugate momentum ir(x), whose evolution 
is governed by the action ( 11.2b . we introduced another set of the canonical variables corresponding to the description in the 
coordinates shifted by a constant dilatation transformation: 

C(x) := C(t, e~ s x) , it(x), (1.7) 

where s is a time independent c-number and n(x) is the conjugate momentum of ((x). In Ref. l27tl . we showed that requesting 
the equivalence between the two quantum systems described by {£, tt} and {£, 7r} guarantees the IR regularity of loop correc- 
tions. Here, the equivalence of two quantum systems means that the same iteration scheme (or formally the same initial condition 
of the interacting system) gives physically the same quantum state in both systems related to each other by the dilatation transfor- 
mation. Namely, all the expectation values evaluated in both systems are equivalent if we take into account how they transform 
under dilatation transformation. Requesting this equivalence will be thought of as the invariance of the initial state under the 
dilatation transformation. In Ref. II2711 . we employed the iteration scheme in which the interaction is turned on at a finite past. 
Then, it turned out that the IR regularity /gauge invariance condition cannot be consistently imposed. In the present paper we 
will set the initial quantum state at the infinite past. We will show that the above transformation can be extended to allow a time 
dependence of the parameter s. As we described in the previous section, this extension plays a crucial role in discussing the 
absence of the secular growth. 

The Euclidean vacuum. The second and third items are related with each other, once we establish the correspondence 
between the gauge invariance and the IR regularity. We will show that the two quantum systems described by {(, tt} and {£, 7f} 
are equivalent if we choose the Euclidean vacuum, which is defined by requesting the regularity of the n-point functions at the 
distant past with the time path rotated toward the complex plane. To be more specific, as the second item, we will show that the 
?i-point functions for ((x) calculated by the canonical variables {£, tt} with the boundary condition of the Euclidean vacuum 
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agrees with the n-point functions for £(f, e s ^x) calculated by the canonical variables {£, 7r} under formally the same boundary 
condition, i.e., 

( C(t, xi)C(t, *a) • • • C(*, *») ) {CM = ( ((t, e'^ Xl )C(t, e'Wx a ) • • • C(*, e >(() a;„) > {c ~ - } . (1.8) 

Combined with the previously mentioned technique to deal with the gauge issue, we will show that when we choose the Eu- 
clidean vacuum, the Hamiltonian density for {£, if}, can be expressed only in terms of the IR irrelevant operators. 

The IR regularity and the absence of the secular growth. As for the third item, we evaluate the n-point function of the 
genuinely gauge invariant operator. Performing the quantization in the canonical system of {£, if}, we will show that the IR 
contributions do not diverge and that the secular growth is suppressed. We carefully investigate the contributions from the 
modes which are comparable to or less than the Hubble scale, i.e., k > aH, without employing the asymptotic expansion with 
respect to k/ aH. As is stressed at the end of the preceding subsection, this point is one of the necessary ingredients to show 
the absence of the secular growth. One may naively expect that the UV modes with k/aH > 1 will not effectively contribute 
to the vertex integration because of the oscillatory behaviour. A more careful consideration tells us that this naive expectation is 
not necessarily correct. In general, vertex integrations become a mixture of the positive and negative frequency mode functions, 
which yields the phase in the UV limit e lr i( k i- k 2+k3 ) wnere r\ represents the conformal time which runs from — oo to 0. Then, 
the phase does not necessarily exhibit the rapid oscillation even for the modes with k m /aH ~ —k m r) > 1, where m = 1, 2, • • ■ , 
which can be a cause of secular growth. Intriguingly, choosing the Euclidean vacuum plays a crucial role not only in the IR limit 
but also in the UV limit. One can show that there is no mixing between the positive and the negative frequency modes, if we 
choose the Euclidean vacuum. Therefore, secular growth is evaded in this case. 

The outline of the paper. The outline of this paper is as follows. In Sec. [II] we will briefly review the way to construct the 
genuinely gauge invariant operator 9 R, following Refs. I25il26ll . Then, we will introduce the canonical variables {(, if} and will 
derive the Hamiltonian for these variables. In Sec. [Hi] we will discuss the items 1 and 2 that we mentioned above. In Sec. 1111 Al 
we will describe the boundary conditions of the Euclidean vacuum and will prove Eq. dl.8t , which implies that the boundary 
conditions of the Euclidean vacuum select the same ground state both in {£, tt} and {£, if}. In Sec. IIII Bl and Sec. MI CI we will 
formulate the canonical quantization in terms of {(, if} and will show that the interacting vertices for these canonical variables 
consist only of the IR irrelevant operators. Particularly in Sec. IIII CI we will show that using the residual gauge degrees of 
freedom, the non-local operator d~ 2 can be made IR regular. In Sec. |IV] we will discuss the item 3. In Sec. II V Al we will 
show that the boundary condition of the Euclidean vacuum leads to the so-called ie prescription in a perturbative expansion. In 
Sec. UV B| we will calculate the Wightman propagator, by which the n-point functions are expanded. Then, in Sec. |IV CI we 
explicitly evaluate n-point functions to investigate the IR regularity and the secular growth. In Sec.[V] as concluding remarks, we 
discuss another possibility of the initial state which satisfies the IR regularity /gauge invariance conditions. We will also mention 
the related papers to clarify what is new in this paper. 

The advantage of the in-in formalism. In our previous publications 1I27I - I29I1 . in calculating n-point functions, we used 
the retarded Green function to solve the non-linear Heisenberg equation. This is because we thought that using the retarded 
Green function, whose Fourier mode is regular in the IR limit, makes the proof of the IR regularity transparent. However, 
the perturbative expansion using the retarded Green function is not suitable for the present purpose, because the positive and 
negative frequency modes are mixed in the vertex integrations once the retarded Green function is used. Therefore, the boundary 
conditions of the Euclidean vacuum does not guarantee the convergence of the time integrations for all the vertices. By contrast, 
when we calculate the n-point functions in the in-in formalism, all vertex integrals can be made manifestly convergent by 
adopting the boundary conditions of the Euclidean vacuum (see Sec. IIII Al ). Since the n-point functions obtained from the 
solution written in terms of the retarded Green function agree with those obtained in the in-in formalism, the vertices which 
do not converge should vanish in the final result of the ?i-point functions. However, the cancellation is obscured in an explicit 
perturbative expansion. Therefore, in this paper, we calculate the n-point function totally based on the in-in formalism, without 
using the retarded Green function. 



II. CONSTRUCTING THE GAUGE INVARIANT QUANTITY 

In this paper, as an explicit model of inflation, we consider a standard single field inflation model whose action takes the form 

M 2 r 

S = 2 / ^ [R ~ ^"^ v ~ 2V ^)] d * x > (2- 1 ) 

where M p \ is the Planck mass and we set to a dimensionless scalar field, dividing it by M p \. However, as long as we consider a 
scalar field with the second-order kinetic term, an extension proceeds in a straightforward way. In Sec. Ill Al we will construct the 
genuine gauge invariant operator corresponding to the spatial curvature of a ^-constant surface. In Sec. Ill Bl we will introduce 
the canonical system {(, if} whose Hamiltonian density is composed only of the IR irrelevant operators. 
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A. Gauge invariant operator and quantization 

We fix the time slicing by adopting the uniform field gauge 6<f> = 0. Under the ADM metric decomposition, which is given by 

ds 2 = -N 2 dt 2 + h ij (dx i + N i dt)(dx j + N j dt) , (2.2) 

we take the spatial metric /iy as 

% = e 2 0>+O[ e *7].. , (23) 
where a := e p is the scale factor, ( is the so-called curvature perturbation and 5%j is a traceless tensor: 

fry*, - . (2.4) 
As spatial gauge conditions we impose the transverse conditions on Sjif 

diSfj = . (2.5) 

Since the time slicing is fixed by the gauge condition 5(f) = 0, there are remaining residual gauge degrees of freedom only in 
choosing the spatial coordinates. In this paper, we neglect the vector and tensor perturbations. The tensor perturbation, which is 
massless, can also contribute to the IR divergence of loop corrections. We will address this issue in our future publication. 

Following Refs. I25II26I1 . we construct a genuine gauge invariant operator, which preserves the gauge invariance in the local 
observable universe. For the construction, we note that the scalar curvature S R, which transforms as a scalar quantity under 
spatial coordinate transformations, becomes genuinely gauge invariant, if we evaluate it in the geodesic normal coordinates on 
each time slice. The geodesic normal coordinates are introduced by solving the spatial three-dimensional geodesic equation: 

d 2 x^ dx^ dx^ 

where a T l jk is the Christoffel symbol with respect to the three dimensional spatial metric on a constant time hypersurface 
and A is the affine parameter. Here we put the index gl on the global coordinates, to reserve the simple notation x for the 
geodesic normal coordinates, which will be mainly used in this paper. We consider the three-dimensional geodesies whose 
affine parameter ranges from A = to 1 with the initial "velocity" given by 



da£,(sc,A) 



(1A 



= e- ax=0) x l . (2.7) 

A=0 



A point x % in the geodesic normal coordinates is identified with the end point of the geodesic, x l g[ (x, A = 1) in the original 
coordinates. Using the geodesic normal coordinates x\ we perturbatively expand x l gl as x gl — x l + 8x l {x). Then, we can 
construct a genuinely gauge invariant variable as 

9 R{t, x) := s R{t, x l gl {x)) = s R{t, x l + 5x l {x))) , (2.8) 

where t denotes the cosmological time. 



B. Dilatation symmetry in the global universe 

The focus of this subsection is on the dilatation transformation, shifting to the rescaled spatial coordinates: 

x l := e s(t) x\ (2.9) 

Solving the Hamiltonian and momentum constraint equations, we can derive the action that is expressed only in terms of the 
curvature perturbation ((x), which is schematically written as 

S = J dtd 3 xL[d t ({x)X(x)}, (2.10) 

Using the curvature perturbation ( and the conjugate momentum defined by it := SL /5{dtQ, the Hamiltonian density is given 
by the Legendre transform as 

5C[C(x),tt(x)] :=Tr(x)d t ((x)-L{d t ((x), ((x)} . (2.11) 
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What is important here is only the fact that the curvature perturbation £ appear in the action either with differentiation or in 
the form of the combination of the physical distance e p+ ^dx l27ll . In the new coordinates ( |2.9t , the physical distance is written 
as e p+< >( t ' xS> ~ s ( t )dx, with the definition of a new variable 



C(t, x) := C(t, x) 



(2.12) 



Thus, if the field ((x) is replaced with ((t, x) —s(t) under the change of the coordinates from x to x, the action basically remains 
invariant. To express dt((x) in terms of the new variable £, we denote the partial differentiation with the spatial coordinates x 
fixed as (d t ((t, x)) x . The subscript associated with the parentheses specifies the spatial coordinates that we fix in taking the 
partial differentiation. Then, we have 



(2.13) 



For brevity, when the fixed spatial coordinates are identical to the ones in the argument of the variable, we simply use dt- Then, 
we can establish an identity 



dtd d xL[d t ax)X(x)} = / dtd d xL[{dtC{t,x)) x X(t,x)^s{t)} 



(2.14) 



Recalling the relation between x and x ( 12.91 ). this equality also means the equality at the level of Lagrangian density, 

e- 3 ^)£[d t C(x), C(x)} = L[(d t at, x)) x , C(t, x) - s(t)}. 

We introduce the canonical conjugate momentum corresponding to ((t, x) in the standard way as 



n(t, a;) 



dL[(d t C(t, x)) x ,C(t, x)-s(t)} 



Noticing the relation 



d(d t c(t, x)) 

d t ({t, x) = (d t ((t, x)) x - s(t) x ■ d x ((t, x) , 



we have 



dL{(d t ((t, x)) x ,C(t, x) -«(*)] 



-3s 



(t) dL[dtt(x),C(x)] _ 3s(t) 



nit, x) = — — '— — — : — - = e "" v "' - = e 

d((d t C(t, x)) x ) d{d t Q{x)) 



ir(x) . 



As is expected, using the commutation relations for £ and ir together with Eqs. ( |2.12t and ( |2.17t , we can verify 

C(i, x), jf(t, y)] = e- 3s Wi5W((x - y)) = i6<V(x - y) , 

as well as 

C(t, i), C(*, V)\ - Ss), 7r(t, y)} = 0. 
The Hamiltonian density for ((x) and tt(x) is obtained in the standard way as 



(2.15) 



(2.16) 



(2.17) 



(2.18) 



(2.19) 



•K 



C(t, x), n(t, x) := %{t, x)d t C(t, x) - £[(d t ((t, x)) x ,((t, x) - s(t)} 

= rt(t, x)(d t ((t, x)) x - L[(dtC(t, x)) x , C(i, x) - s(t)] - s(t) n(t 7 x)x-d x ((t, x) 
= 'KiCit, x) — s(t), Tr(t, x)] — s(t)n(t, x)x ■ d x ((t, x) , 



(2.20) 



where in the equality on the second line we used Eq. (12. 16b . The last equality is exactly the same Legendre transformation as in 
the original system and therefore we can use the same functional form of the Hamiltonian density !K. 

Assuming that s(t) is as small as ((x) and tt(x), we decompose the Hamiltonian densities CK and !K into the non-interacting 
parts, which include only the quadratic terms, and the interacting parts as 



5£[C(z), tt(x)} = J£o[C(ar), n(x)} + 5Cj[C(x), tt(x)} , 



(2.21) 



and 



C{x), Tt(x) 



C(x), tt(x) 



C(x), n(x) 



(2.22) 
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In the above we used the coordinates x instead of x for the {£, 7r} system, but it will not cause any confusion after the relations 
between the {(, tt} and {£, 7r} systems have been established. Here, we replaced 3io[((x) — s(t), tt(x)] with J-Co[((x), tt(x)], 
since ((x) always appears with the spatial derivative in , Ko[((x), tt(x)]. Remarkably, the non-interacting part of the Hamiltonian 
density does not change at all under the dilatation transformation. Using Eq. ( 12.20b , we find that the interaction Hamiltonian 

!Kj[£, tt] is given by 



5£r 



C(x), n(x) 



C{x) - s(t), tt(x) - s(t)ir(x) x-d x C,(x) . 



(2.23) 



In this way, we can write down "Kj only in terms of Q{x) — s(t), ( with differentiation, if and s(t). In Ref. B27I1 . we introduced 
the two sets of the canonical conjugate variables which are connected by the dilatation transformation with a constant parameter 
s. When we take the limit where s(t) is constant, the Hamiltonian density !K(x) takes the same functional form as 3~C(x) 
except for the constant shift of £(x) by — s. It is because, without modifying the gauge condition, we can perform the dilatation 
transformation with the constant parameter s also in the whole universe. Then the action which preserves the diffeomo rphic 
invariance becomes invariant under the change from £(x) to e~ s x) — s. Here we have extended the argument in Ref. Il27ll 
to allow s to depend on time. As we mentioned in Sec. [I] this extension plays the crucial role in our discussion about the secular 
growth. In the next section, we will show that all the interaction vertices in the canonical system {C, tt} are composed only of 
the IR irrelevant operator. 



III. INTERACTION HAMILTONIAN WITH THE IR IRRELEVANT OPERATORS 

In this section, we describe the first two of the three items we raised in Sec. Q] In the preceding section, we derived the 
Hamiltonian for the canonical variables <^(x) and tt(x). Since {C, tt} and {£, tt} are connected by the canonical transformation, 
if we choose the same initial state in both of the two canonical systems, the n-point functions for the same operator, for instance 
9 R, calculated in these canonical systems should agree with each other. However, even if we adopt operationally the same 
scheme to select the initial state in these two systems, it does not guarantee that the selected initial states are the same. In 
Sec. lHI Al after we describe the definition of the Euclidean vacuum, we will show that the condition of the Euclidean vacuum 
operationally selects the same quantum state irrespective of the choice of the canonical variables. This ensures the equivalence 
of these two canonical systems including the choice of the initial quantum state, which we mentioned in the item 2. In Sec. lIII Bl 
we will perform the quantization using the canonical variables {£, tt}. As we will show in Sec. lIII Cl by virtue of the equivalence 
between the two canonical systems, the interaction vertices for {£, 7r} can be expressed in terms of operator products composed 
only of the IR irrelevant operators. 



A. Euclidean vacuum and its uniqueness 

In the case with a massive scalar field in de Sitter spacetime, the boundary condition specified by rotating the time path in 
the complex plane can be understood as requesting the regularity of correlation functions on the Euclidean sphere which can be 
obtained by the analytic continuation from the ones on de Sitter spacetime. The vacuum state thus defined is called Euclidean 
vacuum state. Because of the similarity, here we also refer to the state which is specified by a similar boundary condition as 
the Euclidean vacuum. To be more precise, we define the Euclidean vacuum as follows. In the in-in formalism, the insertion of 
interaction vertices is ordered along the closed time path. By rotating the time path toward the imaginary plane, the forward time 
evolution begins at rj(ti) = — oo(l — ie) and ends at the final time tf and the backward time evolution begins at tf and ends at 
rj(ti) = — oo(l + ie). Here we set e to a small positive number. Since rotating the time path can be better understood by using 
the conformal time n, we introduced the conformal time ?/ as 

We define the Euclidean vacuum, requesting the regularity of the n-point functions with an arbitrary natural number n in the 
limit of rj(ti) — > — oo(l ± ie), i.e., 

F n {xx, ■■■x n ) := (T c C(x 1 )---C(^ n )) <°o as r?(t a ) -oo(l ± ie) , (3.2) 

where a — 1, ■ • • n and T c denotes the time ordering along the closed time path. We first show that the n-point functions of ( 
are uniquely fixed by requesting the condition (13.2b . In this paper, for simplicity, we assume that e p p(p) is rapidly increasing in 
time so that 

\ V (t)\ = 0(l /e^p(t)) . (3.3) 
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Next, we show that the boundary condition of the Euclidean vacuum uniquely determines the n-point functions 
F n (x\, ■ ■ ■ x n ). We schematically describe the Heisenberg equation for £(x) as 

£( = Snl[Q, (3.4) 

where L is the second-order differential operator: 

Z:=d 2 p + (3-e 1 + e 2 )d p --j^. (3.5) 

For notational convenience, we introduced the horizon flow functions, 

£i := — T-7-p, £ n ■= —e n -i, (3.6) 

P dp e n -i dp 

with n > 2, but we do not assume that these functions are small to keep the background evolution unconstrained except for 
requesting Eq. (13.3b , which is valid, for instance, when e„ are constant in time. Using the Heisenberg equation (13 .41 . we can 
obtain the evolution equation of the path-ordered n-point functions F n {x\, ■ ■ ■ x n ) as 

&x a F n {x\, • • • , Xn) = [{-^m}m>7i] 3 (3-7) 

where L Xa is the derivative operator L given in Eq. ( 13.5b with the coordinates x replaced with x a . Since the equation of motion 
for £(x) is non-linear, the equation (13.7b includes the source term (the right hand side) composed of m-point functions of £(x) 
with m > n. We can verify the uniqueness of the n-point functions for £(x) by showing that solution of Eq. ( 13.7b is uniquely 
fixed by the boundary condition (13.2b . To show this uniqueness, we formally solve the equation ( 13.7b as 

F n (xi, ■ ■ ■ , X n ) = /n(xi, ■ ■ ■ , X n ) + &xl*NL [{ F m}m>n] , (3.8) 

where f n (xi, ■ • • , x n ) is a homogeneous solution, while we assume that the specific solution -C^ 1 V^[{F m } m >„] satisfies the 
regularity condition in the limits rj(t a ) — > — oo(l±ie). Now the question is whether the boundary condition ( 13.2b allows us to add 
any homogeneous solutions. In the Fourier space, /„ can be expanded by e~ lfc?? ( ta ) or e %kr)tyta ^ in the limits rj(t a ) — > — oo(l±ie). 
The regularity at r)(t a ) — > — oo(l + it) accepts e~ jfel? ^ ta - 1 only, while the regularity at r](t a ) — > — oo(l — ie) accepts the other. 
Thus the regularity condition in the two limits does not allow to add any homogeneous solutions /„, which implies that the 
n -point functions F n (xi, ■ ■ ■ , x„) are uniquely fixed by the boundary condition of the Euclidean vacuum. 

Next, we show that this uniqueness is ensured independent of whether we use the canonical variables {£, ir} or {£, 7r}. We 
employ the boundary condition of the Euclidean vacuum for the canonical variable £ as well, requesting 

(T c C(xi)---C(x„)) {c>} < oo as ri(t a ) -> -oo(l ± ie) . (3.9) 
Then, we can show that the path-ordered n-point functions 

F„(xi, ■••xn) := (T c C(h, e s ^xi) • • • C(«n, e s ^x n )) { i s} , (3.10) 

agree with the n-point functions F n (xi, ■ ■ ■ x n ) = (T c C(xi) ■ ■ ■ C( x n)){c,7r} fixed by the boundary condition ( 13.2b . i.e., 

F n (xi, ■■ -x n ) = F n (xi, ■ ■ -x n ). (3.11) 

Here putting the suffixes {£, w} or {C, fr}, we denote the canonical variables used in imposing the boundary condition explicitly. 
We again schematically describe the Heisenberg equation for £ as 

L( = S NL [{]. (3.12) 

Since £(x) and ^(x) are connected by the canonical transformation, the equation of motion obtained by operating L on 

C(x) = ((t, e s ^x) = C(x) + s(t)x ■ d x Q{x) + ■ ■ • , (3.13) 

can be recast into Eq. (13.4b by using Eq. d3.12b . A similar argument follows for the equations of motion for the correlation 
functions F n and F n . Using the equation of motion for the n-point functions of ((x), which can be derived from Eq. (13.12b , we 
can confirm that an operation of L Xa on 

F„(xi, ■■■xn) = (T c C(xi)---C(a;r l )) {f)#} + s(ii)(T c xi-5 a!1 C(ti, xi)--.C(*», *n)} { ^ } + • • ■ (3.14) 



10 



leads to 

£ Xa F n (xi, ■■■ , x n ) = V$ L [{F m } m>n ] . (3.15) 

This equation takes the same form as the equation of motion (13. 7b . We also note that the boundary condition of the Euclidean 
vacuum (13.91 ) implies 

F n (xi, ■■■ , x a ,-- ■ , x„) < oo as - r](t a ) ->• oo(l ± ie) . (3.16) 

The equivalence ( 13.111 ) is now transparent, because the equations of motion ( 13.71 ) and ( 13.151 ). and the boundary conditions ( 13.21 ) 
and ( 13.161 ) are the same, and the latter specify the solutions of the former uniquely. This equivalence is a distinctive property of 
the Euclidean vacuum 1 . Here we took the boundary conditions for n-point functions as the definition of the Euclidean vacuum 
state, assuming the existence of such a quantum state. In Sec. lIV Al we explain such a Euclidean vacuum, if exists, should be 
the one given by the ordinary ie prescription. 



B. Rewriting the n-point functions 

In this subsection, we rearrange the expression for the n-point functions of the genuinely gauge invariant variable 9 R into a 
more suitable form to examine the regularity of the IR contributions. First, solving the three dimensional geodesic equations, we 
obtain the relation between the global coordinates x % gl and the geodesic normal coordinates x 1 as 



<</' 



-C(t,e Cm) x i + ... f (3 1?) 



where the ellipsis means the terms which vanish when £(x) is spatially homogeneous, i.e., the terms suppressed in the IR limit. 
Note that changing the spatial coordinates into the geodesic normal coordinates also modifies the UV contributions. Tsamis and 
Woodard l45ll showed that using the geodesic normal coordinates can introduce an additional origin of UV divergence, which 
may not be able to be renormalized by local counter terms ll46ll . It should be clarified whether this issue is a serious problem or 
not, but we defer it to a future study. Instead, to keep the UV contributions under control, we replace C(x) in Eq. (13. 17b with the 
smeared curvature perturbation 9 ((t), i.e., 



l 9l 



sf(t V , (3.18) 



with 



%(t) := J - r A3 u , " , (3-19) 

J d 3 xW Lt {x) 

where Wh t (x) is a window function which is non-vanishing only in the local region S t n 0. We approximate the averaging 
scale at each time t by the Hubble scale, i.e., L t ~ 1 / {e pi ^ p(t)} . Although 9 Q appears on the right-hand side of Eq. (13.19) , 9 ( 
is defined iteratively at each order of the perturbation. We calculate the n-point functions of 3l x 9 ((t, x), instead of 9 R, with 

9 C(t, x) := C(t, e- sm x) . (3.20) 
Here, "R x denotes the IR suppressing operator such as 

d x f J d 3 xW Lt (x) \ 

e^wy y 1 - sd* y w Lt (v) r '"- (3 - 21) 

where x is the spacetime coordinates of the field on which these operators act. Although 3l x 9 ((t, x) is not genuinely gauge 
invariant, it is still invariant under the dilatation transformation, which is associated with the dominant IR contributions. In fact, 
since the smeared curvature perturbation 9 ((t) transforms into %(t) — f under the dilatation transformation: x — > e~?x with a 
constant /, ^^(x) is kept invariant under this transformation. By contrast, the constant part of 9 ((x) can be modified under the 
dilatation transformation as 9 ((x) — > 9 C(^) — /■ Since the genuine gauge invariant variable 9 R(x) should not be affected by the 
dilatation transformation, which is a part of the residual gauge transformations, %(x) appears only in the form of 3l x 9 ((x) when 



The uniqueness of the Euclidean vacuum becomes intuitively clear when the Hamiltonian is time independent and the lowest energy eigenstate is non- 
degenerate, because the ie prescription selects the unique ground state of the system. 
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we express g R(x) in terms of 9 ((x). As we can compute 9 R(x) from % x 9 ((x), our goal is to prove that the expectation values of 
products of Jl x 9 £(x) are IR regular. 

First, we calculate the n-point functions of 9 Q without the IR suppressing operator 3l x : 



(0\%(t f ,x 1 )---%(t f ,x n )\0). 
Using the eigenstates of 9 ((tf), \ s) H which satisfy 9 ((tf) \ s) H — s \ s) H , we can construct a unit operator 



ds | s) H H (s | 



(3.22) 



(3.23) 



Inserting it into the expression for the n-point functions, we obtain 

(0\%(t f , *!)••• %(t f , x n )\0) = J ds( 0|C(*/, e-'xi) ■■<(</, e- s x n )\ s) Js\o ) 



ds( 



C(tf, xi)--C(t/, x n ) 



S) ( s 

H H 



(3.24) 



In the first line we could simply replace 9 C(tf) with s, because 9 ((tf) and C(t/, a; ) commute with each other. Since the 
Heisenberg picture field x) is related to the interaction picture field (i(t, x) as 

C(t, x) = U\{t)Q{t, x)U T (t), 



(3.25) 



where the unitary operator Uj (t) is given by 



U^t) := lim Texp 

r)(ti)— > — oo(l— ie) 



rj (t i ) — > — oo ( 1 — ie) 



n=0 



dt I d 3 x& T 



t-n. 



(l(x), 7T/(x) 



lim 5Z(-i) n / dt„ / dt n _i--- / dii 



x / d 3 a;„ • • • / d 3 a3i 5£j 



Cj(ie„), ffj(a;„) 



■ 5Cj 



(3.26) 



Thus, the n-point function can be rewritten as 



(0\%(t f , Xl )---%(t f , x n )\0) = ds(0 



Texp 



i / dtd^fti 



x (i(t f , xi)-- -C/(*/) a: n )Texp 



C(x), 7r(x) 

-i / did^IK/ 



C(x), 7r(x) 



H H 



(3.27) 



where T denotes the anti time-ordered product. Notice that the interaction Hamiltonian JCj does not contain the second or higher 
derivative of s(t). We construct unit operators, using the eigenstates |s(i)) and \s(t)) which satisfy 



where 



9 ut)\s(t) > = *(*)!*(*)>, 3 ut)\m > = m\m 

J d*xW Li {x)Q{t,e- a ZiMx) 



Jd 3 xW Lt (x) 



(3.28) 



(3.29) 



is the smeared interaction picture field. We next replace all s(t) and s(t) with 9 Ci(t) and 9 (j(t), respectively, by inserting the 
unit operators; 



ds(t) s(t))( S (t) 



dm m)(m 



(3.30) 

To perform this replacement without ambiguity, we fix the operator ordering in JCj to the Weyl ordering, in which (i(x)—s(t) and 
tti(x) are symmetrized. Instead of considering the explicit form of the interaction Hamiltonian, we use a schematic expression 
of 5ii which is expanded in a power series of s(t) as 



Cr(x), 7r/(i)l - %i \ti(x) - s(t), 7ri(x)] = J2 \Ci(x) ~ s(t), ni(x) 



a=0 



(3.31) 
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although a is at most 1. Here, we stress that the perturbations (i(x) and s(t) appear in the Hamiltonian density "K[ only in the 
form of (i(x) — s(t) or its spatial differentiations. Inserting the unit operators, we obtain 



Cr(aO, 7fj(a:)] = E / ds W / ds(t) {s(t)} a ^C (a) [&(*)-*(*), 



(3.32) 



After we replace s(i) with 9 (j(t), 9 (j(t) is located next to the operator \s(t))(s(t)\. Noticing the fact that s(t)\ s(t) ) can be 
expressed as 



s(t)\s(t))=% I (t)\s(t)) = J d 3 xW Lt (x) 1 ( ))j 

where in the second equality, we replaced 9 (j(t) in the argument of Q with s(t), we use 9 (j(t) expressed as 

g - fd 3 xW Lt (x)<ii{x) 
^ I[> Jd 3 xW Lt (x) ' 

instead of the expression given in Eq. (13.291 l, when we replace s(t) with 9 (j(t). Using the formula 

%(t, x) - a(tj) A\ s(t) ) = (Cr(t, x) - %j(t)) A\ s(t) ) + [%j(t), A] \ s(t) ) , 



(3.33) 



(3.34) 



(3.35) 



we replace (d(x) — s(t)) with (Ci(x) — 9 (i(t)) one by one. By induction, the operator A is supposed to be composed of 
(i(x) — s(t) and tti(x). Since 9 (j(t) commute with (i(x) — s(t), the non-vanishing commutation relation is only the following: 



[<&(*),*/(<,*)] 



Jd 3 x W Lt (x) 



d 3 yW Lt {y) 0(t, x), 7f/(t, y) 



'/d 3 a! W Lt (a!) ' 



where we used 



Ci{t,x), 7r/(t, y)l = C7}(t) [C(t,x), 7r(t, y)l ^(t) = «( 8 )(x - y) . 



(3.36) 



(3.37) 



Since the commutator including 9 Ci(t) yields only a local function, we can conclude that operators left after exchanging s(t) 
with 9 Ci(t) are also composed of (i(x) — s(t) and tti(x). Repeating this procedure, we can replace all s(t) with 9 (j(t) as 



ct=0 



(3.38) 



where to denote the modification after the replacement of s(t) with 9 (j(t), we put ' on the interaction Hamiltonian. Replacing 
s(t) with 9 £(t), we obtain 



C/(x), nj(x)] = E / ds W / [Ci(x)-%At), 5f z (x)] |«(t))(s(t)|*(t))(«(t)| { 9 C/W}" • (3-39) 



We repeat this procedure for all integrating Hamiltonian densities which appear in the perturbative expansion of the n-point 
functions ( 13.27b . After these replacements, the possible dependence of the n-point functions on s(t) and s(t) remains only in 
\s(t))(s(t)\ and \s(t))(s(t)\. Since requesting the Euclidean vacuum uniquely determines the initial state independent of s(t) 
and s(t), we can remove the identity operators J ds(t) \s(t))(s(t)\ and J ds(t) \s(t))(s(t)\ as long as we choose the Euclidean 
vacuum. (From the same argument, we can remove the identity operator J ds \ s) H H (s\.) Then, the Hamiltonian density is recast 
into 



3^J 



Ci(x), ttjCx)] ^ E^W [Ci(x)~ 9 Ci(t), nj(x)] { 9 Ut)Y ■ 



(3.40) 



a=0 



Note that we can express 9 (j(t) as 



9 Ut) = J d 3 xd t < 


\ W Lt (x) \ 


Jd 3 xW Lt {x) J 


= J d 3 x d t < 


\ W Lt (x) \ 


Jd 3 xW Lt (x) } 



Ci(x) + 



J d 3 xW Lt (x)d t Ci(x) 
J d 3 xW Lt (x) 

} Jd 3 xW Lt (x)d t (j(x) 
Jd 3 xW Lt (x) 



(3.41) 
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where in the last equality, we inserted = 9 d(t) dt{J d 3 xWL t (x )/ J d 3 xWi, t (x)} and the last term in the last line can be 
written in terms of ttj(x). 

In this way, we can show that all (js in the interaction vertices are multiplied by an IR suppressing operator H x . Notice that, 
replacing the c-number parameter s(t) with the operator 9 (j(t), we rewrote the Hamiltonian density as in Eq. (13.401 . In this 
procedure, we used the fact that the initial state specified by the boundary condition of the Euclidean vacuum does not depend 
on the choice of the canonical variables. We should emphasize that if this equivalence of the initial state were not guaranteed, 
we could not express the interaction Hamiltonian only in terms of (js with an IR suppressing operator. 

C. Restricting the interaction vertices to the local region 

In the above discussion, we found that the interaction picture fields which appear in the interaction vertices can be expressed 
only in terms of ttj(x) and Ci( x ) ~ 9 Ci(t)- Now, we can verify the item 1 presented in Sec. [I] which claims that the interaction 
vertices are constructed only from the IR irrelevant operators. As we showed in the previous subsection, all the interaction picture 
fields are associated with an IR suppressing operator Ck x , which increases the power law index with respect to the wavenumber 
k in the IR limit. To complete the proof of the argument given in the item 1, we need to show that the inverse Laplacian d~ 2 , 
which appears in solving the constraint equations to obtain the lapse function and the shift vector, does not reduce the power 
law index with respect to k in the IR limit. The potential danger can be understood as follows. When we choose the boundary 
condition specified by the regularity at the spatial infinity following the standard procedure, the action of the operator d~ 2 yields 
a multiplicative factor 1 /k 2 . This IR singular behavior arises because the information from the outside of our observable region 
is used to determine the lapse function and the shift vector. 

To remove this potential IR singular behavior originating from the inverse Laplacian, we need to discuss the causality. The 
causality is basically maintained even at the quantum level in the sense that the interaction vertices located outside our observable 
region are decoupled in the in-in formalism. In the ordinary field theory with a local interaction, this can be shown by system- 
atically replacing the Wightman function G + with the retarded Green function plus G~ (see Appendix of Ref. 12711 ). However, 
when the gravitational perturbation is taken into account, it becomes less transparent whether the causality is maintained owing 
to the issue of the lapse function and the shift vector mentioned above. 

Here, we should recall that what we really need to evaluate is the expectation values of genuinely gauge invariant variables, 
which do not depend on the choice of the residual gauge degrees of freedom. As we explicitly showed in Appendix lA] using the 
residual gauge degrees of freedom, we can modify the boundary conditions of the lapse function A*" and the shift vector N so 
that the terms associated with d~ 2 are completely specified by the fields within the local region 0. Then, the operation of the 
non-local operator d~ 2 no longer reduces the power law index with respect to k. In this way, using the degrees of freedom in the 
choice of boundary conditions, we can localize all the interaction vertices within the causally connected local region 0. Since 
$-x((x) is not invariant under the residual gauge transformations, their ?i-point functions are not invariant in general under the 
change of boundary conditions of N and Ni. However, when we calculate n-point functions for the genuinely gauge invariant 
operator 9 R using those for 3l x 3 C> changing the boundary conditions should not affect the result. 

IV. THE IR REGULARITY AND THE ABSENCE OF THE SECULAR GROWTH 

In this section, we will calculate the n-point functions of 3i x 9 ^(tf, x), properly taking into account not only the IR-modes 
but also the modes with k\rj(t) | > 1. As stressed in Sec. [I] to prove the absence of the secular growth, we need to evaluate the 
contribution of the latter modes carefully. In the preceding section, we showed that, using the canonical variables and w(x), 
we can expand the n-point functions of H x 9 ((x) for the Euclidean vacuum only in terms of the IR irrelevant operators. In this 
section, based on the perturbative expansion in the {(, it} system, we will discuss the IR regularity and the absence of secular 
growth in the n-point functions. 

For our current discussion, the explicit form of the interaction Hamiltonian density !Kj is not necessary. We use a formal 
expression 

n 

Kj[Cr(x), 7f/(x)] = M p Wei (*) E A W IT ^i m) Ci( x ) , (4.1) 

n— 3 m—1 

where A(t) is an 0(1) dimensionless time dependent function which can be expressed only in terms of the horizon flow functions. 
To discriminate different IR suppressing operators, we associate a superscript (m) on H x . 
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A. Euclidean vacuum as is obtained by ie prescription 



In the preceding section, we introduced the Euclidean vacuum as a vacuum state which satisfies the boundary condition 
(I3.2t / (|3.9l l. Here we show that this condition forces us to adopt the ordinary perturbative description of the ie prescription. We 
expand the curvature perturbation Q (x) as 



Cz(s) 



where a k is the annihilation operator, which satisfies 



d 3 fc 

(2tt) 3 / 2 



lkx v k {t)h k +h.c. 



a k , al, 



= S^(k-k'), [a k ,a k }=0. 



The mode function v k (t) should satisfy 



dl + (3 - e x + e 2 )d p 



k 
ePp 



v k = 



(4.2) 



(4.3) 



(4.4) 



Since the boundary condition (I3.2b / (|3.9b should hold at the tree level, the asymptotic form of the positive frequency mode 
function v k (t) should be oc e^ tkn ^> . Factoring out this time dependence at r\ — >• — oo, we express v k (t) as 



v k (t) 



At) 

£3/2 



where we introduced 



A(t) 



fk(t)e 



P(t) 



-ikr](t) 



Vei(t)M p i 



(4.5) 



(4.6) 



as an approximate amplitude of the fluctuation. The function f k (t) satisfies the regular second order differential equation with 
the boundary condition 



fk(t) -> 



V2ePp 



for 



k V (t) 



(4.7) 



Since both the differential equation and the boundary condition of f k (t) are analytic in k for any t, the resulting function should 
be analytic as well. Namely, f k (t) does not have any singularity such as a pole on the complex fc-plane. We suppose that a 
positive frequency function for a general vacuum except for the Euclidean vacuum is given by a linear combination of v k and 
vl with the Bogoliubov coefficients which have some nontrivial structure of singularities in the complex fc-plane or diverge at 
infinity. The only exception to evade the singularity is setting the Bogoliubov coefficients to constants, but then the UV behavior 
does not agree with the one in the Minkowski vacuum. 

On the other hand, in the limit —kr](t k ) <C 1, the function f k (t) is proportional to A(t k )/A(t), where t k is the Hubble 
crossing time defined by — kr)(t k ) = 1, because the curvature perturbation should be constant in this limit. Hence, the expansion 
for small fc is in general given by 



A(t)f k (t)=A(t k )[l + 0(k\ri(t)\)] . 
By using Eq. ( |4.5t , the Wightman function is given by 

d 3 fc 



(4.8) 



G + (x, x') 



D ik-(x — x') 



(2tt) 3 
A(t)A(t') 



v k (t)v* k (t' 



d 3 fc 1 

(2tt) 3 fc 3 ' 



fk(m(t')e 



*U'\JHv(t')-r)(t)) 



(4.9) 



Using the in-in formalism, the n-point functions can be expanded by the Wightman function. At this point, the vertex integrals 
should start with r\ = — oo to be able to impose the boundary condition of the Euclidean vacuum d3.2b /( l3T9l >. Although the 
integrands of the vertex integrals are infinitely oscillating in the limit rj —> — oo, the time integration can be made convergent by 
adding a small imaginary part to the time coordinate, which is nothing but the ordinary ie prescription. To see the convergence 
of the time integration more explicitly, we first consider the integral for the vertex which is closest to the past infinity r) 
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v<»{t: {x m }) vw(t';{x m }, { Xm }) 




FIG. 1 : We perform the vertex integrations from the closer vertices to the future or past end of the closed time pass. The left figure represents 
the integration about the vertex which is closest to the past end and the right figure represents the integration about the next to the closest to 
the past end. 



—oo(l — ie) (see Fig. Q]). The interaction picture fields (i(x) included in this vertex are contracted with (i(x m ) contained 
in vertices labelled by m = 1, 2, • • ■ , n, and give the Wightman function G + (x m , x). Then, the vertex integration with n 
interaction picture fields is given by 

V {1 \t', {x m }) := M 2 pl f dt f d 3 xe 3 "W ei (t)p(t) a A(t) f[ X Xm %W G + (x m , x) . (4.10) 

J ** J m=l 

The Euclidean vacuum condition (!3.2b /( f3T9T > requires the convergence of this integral when we send r)(U) —> —00. Since the 
Wightman functions contain the exponential factor e M? w5D"> fcm , the integral can be made convergent by changing the integration 
contour as shown in the left panel of Fig. Q] which is exactly what is known as the ie prescription. 
The vertex integration next to the closest to the past infinity 

VW(t",{x m },{x m ,}) 

dt' I dVe 3 ^'> ei (i')/K0 2 A(i') II ^ m M x m,) G + (x m ,,x')V^(t', {x m }), (4.11) 

' i ra' — l 

can be done in a similar manner, where n' is the number of propagators connecting between this second vertex and the vertices 
other than the first one. If we assume the integration over the time coordinate of the first vertex t up to t', the exponential factor 
in G + [x m , x) can be replaced as 

e ik m (n(t)-ri(t m )) _^ gik OT (»K*')-*K*m)) _ (4.12) 

Therefore all the Wightman functions connecting the vertices at t' or before t' with the vertices after t' give an exponential factor 
which is suppressed by adding +ie to rj. This is again consistent with the boundary condition of the Euclidean vacuum. The 
same argument can be made for the other vertices as well. 



B. The IR/UV suppressed Wightman function 

Since all (i(x)s in the interaction Hamiltonian are multiplied by the IR suppressing operators 5l x , the n-point function of 
Ji x %(x) can be expanded by the Wightman function 'R x 'R x 'G + (x, x') and its complex conjugate H x 3l x iG~ (x, x'). In this 
subsection, we calculate the Wightman functions multiplied by the IR suppressing operator, H x 3l x i G + (x , x') for t > t'. After 
integration over the angular part of the momentum, the Wightman function 3i x Ji x i G + (x, x') can be expressed as 



1 Z" 00 dA- 

O^K^G+ix, x>) = — — % x 'R x> A{t)f k {t)A{t')r k {t') 
In 1 J k 



ik(a+(x, x') — a- (x, x')) 



(4.13) 



where we introduced 



a±(x, x') := r){t') - V (t) ± |: 



We first show the regularity of the k integration in Eq. (14.13b . Since the function fkif) is not singular, the regularity can be 
verified if the integration converges both in the IR and UV limits. The regularity in the IR limit is guaranteed by the presence 
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of the IR suppressing operator. The IR suppressing operators 3i x add at least one extra factor of k\rj(t) | or eliminate the leading 
i-independent term in the IR limit, and yield 



R x A(t)f k (t) 



gika + (x ,x' ) ^iko-{x,x ) 

£fc(<7+(ar, x') — a- (x, x')) 



= A(t k )e ik ^O(k\r,(t)\) 

= A(t)e ikr >^ f{fcb7(t)|} ( " a+1)/2 



(4.14) 



where we have introduced the spectral index n s — 1 := dlog(|.A(ffc)| 2 )/dlogfc. Thus, the operation of "R x makes the k 
integration in Eq. (14.131 regular in the IR limit. Next, we consider the convergence in the UV limit. In Eq. (14.131 . the integration 
contour of k should be appropriately modified at k — > oo so that the integral becomes convergent. This modification of the 
integration contour can be also understood as a part of the ie prescription, because adding a small imaginary part to all the time 
coordinates as n — s- n x (1 — ie) leads to the replacement rj(t') — 77(f) — > r)(t') — 77(f) + ie, where we note n(f') — 77(f) < 0, 
and hence to introducing an exponential suppression factor for large k. This UV regulator makes the integral finite for the large 
k contribution except for the case a± (x, x') = 0, where x and x' are mutually light-like. Since the expression of the Wightman 
function obtained after the k integration is independent of the value of e, the regulator makes the UV contributions convergent 
even after e is sent to zero. For a±(x, x') = 0, the integral becomes divergent in the limit e — > 0, but the divergence related to 
the behavior of the Wightman functions in this limit is to be interpreted as the ordinary UV divergences, whose contribution to 
the vertex integrals must be renormalized by introducing local counter terms. Thus, the Wightman function 3i x 'Ji X 'G ± (x, x') is 
now shown to be a regular function. 

Since the amplitude of the Wightman function with the IR suppressing operator is bounded from above, we can show the 
regularity of the n-point functions, if the non-vanishing support of the integrands of the vertex integrals is effectively restricted 
to a finite spacetime region. Since the causality has been established with the aid of the residual gauge degrees of freedom, 
the question to address is whether vertexes at the distant past is shut off or not. To address the presence of such a long-term 
correlation, we discuss the asymptotic behavior of the Wightman function 'R X 'R X ' G ± (x, x'), sending if to a distant past. Recall 
that when a±(x,x') 7^ 0, we can rotate the integration contour with respect to k even toward the direction parallel to the 
imaginary axis. Rotating the direction of the path appropriately depending on the sign of a±(x, x'), the integrand becomes 
an exponentially decaying function of k. This rotation of the integration contour can be done without hitting any singularity 
in the complex /c -plane, because the function f k (t) is guaranteed to be analytic by construction. If we choose other vacua, 
this operation induces extra contributions from singularities. Since we send f' to the past infinity, assuming \r)(t')\ |n(£)|, 
cr±(x, x') is 0(|n(f')|), except for the region where the two points are mutually light-like 2 . Then, the integration of k on the 
right-hand side of Eq. ( 14- .13b is totally dominated by the wavenumbers with k < l/|j7(f')| -C 1/ 1 77(f) |. Using Eq. ( |4.14t which 
gives the asymptotic expansion in the limit k\rj(t)\ <C 1, we obtain 



R x X. x ,G + (x, x') =A(t)0 



dk 



o k {ePWp(t) 



(n. + l)/2 



% x ,A(t')t k {t'y k ^ 



WW* (ML)' 



(4.15) 



where on the second equality, we performed the k integration, rotating the integration contour. We should emphasize that we did 
not employ the long wavelength approximation regarding the Hubble scale at f' to properly evaluate the modes k of 0(l/|?/(t')|) 
as well. 



C. The secular growth 

In this subsection, focusing on the long-term correlation, we discuss the convergence of the vertex integrals of the n-point 
functions for the Euclidean vacuum. We start with the integration of the n-point interaction vertex which is the closest to 



2 Let's introduce a physical length scale Atj v to remove the contributions from the vicinity of the lightcone. On the time slice specified by 77' , we neglect the 
region within the distance Atjv from the intersection of the light cone emanating from x with this time slice. Under this restriction, we have + (x, x') > 
|?7(t')| i?(t')Auv an d hence rr^lx, x') turns out to grow in proportion to |r?(t')|. This argument might be too heuristic, but we believe that the contribution 
from the region neglected here will not change our discussion about the IR regularity of the n-point functions. In order to clarify this point, it would be 
necessary to incorporate the discussion about the UV renormalization, which is beyond the scope of this paper. 
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IR suppression 



ln(e'p) 



ln(e"o) = lni 




(a) M~l 



UV suppression 



IR suppression 




\n(e"p) = ]nk 



\n(e"p) = hik+ lne 
^Ink 



(b) A/ »/ 



UV suppression 



FIG. 2: These figures show which modes can contribute to the loop integrals in the n-point function of % for the Euclidean vacuum. The 
horizontal axis represents the wavenumber lnfc and the vertical axis represents the time ln(e p p) ~ ln(l/|»j|), which becomes the number of 
e-folding in the limit e\ <C 1. The red region is suppressed because of the operation of the IR suppressing operator % x and the blue region 
is suppressed because of the exponential suppression of the it prescription. The dotted line with log(e p p) = log k is the mode of the Hubble 
scale. The left figure (a) is for the case with M ~ 1 and the right figure (b) is for the case with M ^> 1 . 



r/ = — oo(l — ie). By inserting the expression of the Wightman function 'R x '3l x i G + (x , x') with t ^s> t\ given in Eq. ( 14.151 ) into 
Eq. ( 14. lOt . the vertex integral can be estimated as 



V (1 \t\ {x m }) = 



Ml f At f ^xe^e^pitfmimr f[ MUn) ( 

Jt . J m=l 



(4.16) 



As we have explained in Sec. IIII Cl the interaction vertices are confined within the observable region, i.e., the non- vanishing 
support of the integrand is bounded by \x\ < L t ~ \v(t)\- Thus, we obtain 



VP(t\ {x m }) = 



(4.17) 



As we have performed momentum integral first, the exponential suppression for large \r/\ is not remaining any more. However, 
picking up //-dependence of the integrand of Eq. ( |4.17t , we still find that the contribution from the distant past is suppressed if 



— > as rj — > — oo 



(4.18) 



Then, the time integral converges, and the amplitude of Vn (v'i { x m}) is estimated by the value of the integrand at the upper 
end of the integration as 



vW(t', {x m }) = 



X(t'){A(t')} n ' 2 l[ A(t r 



m—1 



v(t m ) 



(4.19) 



Therefore, when a Wightman propagator is connected to a vertex located in the future of x', i.e., when t m > t', the t-integration 
yields the suppression factor {rj(t m ) /^(t')} 11 ^— . We denote the number of such propagators by h. 
Similarly, we can evaluate the amplitude of as 



V^(t", {x m },{x m ,}) 
= 



-±X'{ v '){A{ v ')r'- 2 ] [ A(tr, 

r l Jt, 



Extracting the //-dependent part in the above expression, we obtain 

rv(t") dr] , 



v(t m ) 

rf 



V^(t(r,'),{x m }) 



! r XW)X'(r ] '){A(r 1 ')r +n ^W\ 



-(n'+h) 



(4.20) 



(4.21) 
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Notice that all the Wightman propagators which are connected to the field Q located in the future of x' yield the suppression 
factor \r}{t")\~ J1J ^ ± . 

Now the generalization becomes easy. For the N v -lh vertex, the temporal integration becomes 

^A^J^wJ^-^lwJ-^', (4.22) 

where Nf denotes the number of £/s contained in the vertices up to the A^-th, M denotes the number of the Wightman propa- 
gators connected to a vertex with r\ > r)N v , and A denotes the product of the interaction coefficient up to the A^-th vertex. Thus, 
the convergence condition is given by 



N f -2N 



as r) -oo . (4.23) 



Since all interaction vertices have at least one Wightman propagator connected with their future vertices, M should satisfy 
M > 1. 

As a simple example, we consider the case where E\ is constant. In this case, A is expressed only in terms of e\ and takes a 
constant value. By assuming M = 1 and using n s — 1 = — 2ei, the convergence condition yields 

-eiJV + (1 - £i) 2 > , (4.24) 

with A" := Nf — 2N V . In the slow roll limit ei -C 1, the above condition is recast into 

A"<0(l/ei). (4.25) 

The intuitive understanding of the above suppression mechanism is as follows. In the Euclidean vacuum case, only the 
contributions around the Hubble scale at each time are left unsuppressed (as shown in Fig. |2j. When only the modes around the 
Hubble scale, i.e., k\rj\ ~ k/e p p = 0(1), are relevant, the Wightman function % x 3l x r G + (x, x') is necessarily suppressed when 
v(t)/v(t') ^ 1- This is because if x and x' are largely separated in time, any Fourier mode in the Wightman function cannot 
be of order of the Hubble scale simultaneously at t and t'. When we consider the contribution of vertices located far in the 
past, at least one Wightman function should satisfy r)(t)/r)(t') <C 1, and therefore it is suppressed. However, when we consider 
a diagram for which a cluster of vertices in a distant past is connected to the vertices around the observation time by a single 
propagator, i.e., in the case with M = 1, the IR suppression comes only from this propagator. When the number of operators 
in the cluster of vertices in the past is sufficiently large, the suppression due to this propagator can be overwhelmed by the large 
amplitude of the fluctuation, which increases as the energy scale of inflation increases in the past direction. This corresponds 
to the case when the condition ( 14.231 ) is broken. However, we should also stress that the contributions from the distant past are 
suppressed and the secular growth never appear in the slow roll inflation, unless the order of perturbative expansion N takes an 
extremely large value such as 1/ei ~ O(10 2 ). When the convergence condition ( I4.23l l is satisfied, all the time integrations are 
dominated by the contributions near its upper end. The order of magnitude of the rc-point functions of 3? X 9 C(^ fi x)is then given 
by 

(Ol^^t/, Xl )R X2 %(tf, x 2 )---% Xn %(tf, x n )\0) ~ X(tf){A(tf)} N . (4.26) 



V. CONCLUSION AND DISCUSSION 



A. Euclidean vacuum satisfies the strong constraint on the initial states 



In this paper, we showed that when we choose the Euclidean vacuum as the initial state, the vertex integration in the n- 
point functions for the genuinely gauge invariant curvature perturbation is regular unless a very high order in the perturbative 
expansion is concerned. Figure|3]shows the outline of the proof. We should emphasize that the regularity of the ?i-point functions 
in the limits 77 — > — 00 (1 ± ie) plays a crucial role in the proof: (i) Requesting this regularity guarantees the equivalence between 
two quantum systems, i.e., the original system in which the Hamiltonian contains the IR relevant operators and the quantum 
system in which the Hamiltonian is totally composed of IR irrelevant operators, (ii) It guarantees the analyticity of the mode 
function Vk (t) with respect to the wavenumber k for arbitrary t. By virtue of the aspect (i), we can rewrite the n-point functions 
of % into those expressed in {£, tt}, in which all the field operators are manifestly associated with the IR suppressing operators, 
3l x . The aspect (ii) leads to the exponential suppression in the UV so that the non-vanishing support of the ^--integration is 
restricted to —ki] < 0(1). It might be intriguing that choosing the Euclidean vacuum plays the crucial role in discussing the 
suppressions both in the IR and UV components. Since these suppressions make the Wightman function (in the position space) 
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Choosing the Euclidean vacuum 
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Regularity of w-point functions for Q and g 
in the limits r\ — >-oo(l+ie) 
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Sec.IIIA, IVA 



(i) Equivalence between the 2 canonical systems 

(ii) Analyticity of the mode function 

Sec.IVB 



Convergence of ^-integration in Wightman functions 
(i) — >IR suppression, (ii) — >UV suppression 



Preserving the gauge invariance 
in the local universe 



Sec.III C 



Interaction verteices are localized 
within the past lightcone. 
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\ 

Sec.V S 



Sec.IVC 

Suppression of the long-term correlation 



— Regularity of vertex integration in the «-point functions 



FIG. 3: The outline of the proof which shows the regularity of the n-point functions of the genuinely gauge invariant variable for the Euclidean 
vacuum. Since we have left a possibility that the n-point functions can become regular without requesting the boundary condition of the 
Euclidean vacuum, we used the dotted arrow. 



associated with an IR suppressing operator regular everywhere except for the light cone limit, the missing piece to prove the 
regularity of the n-point functions is to show that the integration region of each vertex integral is effectively confined to a finite 
portion of the spacetime. Using the residual gauge degrees of freedom, we can confine the interaction vertices within the past 
light cone. Since the long-term correlation is shut off because of the suppression both in the IR and UV, the integration region of 
the vertex integrals is ensured to be effectively finite. Therefore, the n-point functions for the Euclidean vacuum are expressed 
by integrals whose integrand and integration region are both finite, and hence they are manifestly regular. Thus, we conclude 
that the Euclidean vacuum is a suitable initial state of the universe which is free from the IR pathology even in the presence of 
non-linear interactions. 

In this section, we further address the converse question; "When we request that the n-point functions are finite and free from 
the secular growth, is the Euclidean vacuum the unique possible initial quantum state?" To be precise, the condition we impose 
here is the regularity of n -point functions on the real time axis including the distant past, i.e., —77 ^> 1. We naively expect that 
in this case, the Euclidean vacuum is the unique possibility. Since any excitations are blue shifted at an earlier time, any small 
deviation from the Euclidean vacuum state at a finite time will lead to some singular behavior in the limit — 77 ^> 1. However, 
we do not have any rigorous proof about this argument yet. There might be a fundamental obstacle when we try to make this 
statement precise. When we trace back the history of the universe, it should inevitably enter the regime in which the background 
energy density and hence the amplitude of the vacuum fluctuation are so high that the perturbative analysis would not make 
sense any more. 

As an alternative setup of the problem is to require the regularity of the n-point functions just for 77 > 77, with a certain initial 
time, rji. The relaxed requirement of the regularity allows us to take other states, if correlation functions for these states can be 
reinterpreted as correlation functions for the Euclidean vacuum. We introduce a new operator 

A\ m) = J d'x^^OW, 

with an arbitrary choice of the IR suppression operator 3?( m ) where m is just the label for distinction. Then, we can define the 1 
particle state by | 1( TO )) := NA'^ |0) with an appropriate normalization factor N. The n-point functions of 3i x 9 ((x) at the initial 
time ri — rji for the 1 particle state | l( m )) defined at the initial time can be expressed in terms of the (n + 2)-point functions 
for the products of Ji x 9 ((x) for the Euclidean vacuum. When the initial distribution is regular, as we showed in Ref. ll28ll . the 
distribution at late times will be kept regular as well. Similarly, we can construct excited states with plural particles. (Similar 
excited states are discussed in de Sitter spacetime in Ref. JiU].) However, here the allowed number of inserted operators might 
be bounded because our proof of regularity does not apply when the order of perturbation becomes very high. 
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To extend the above discussion to the n-point functions at a later time, we only need to show the regularity of the n-point 
functions which are defined as the expectation values of the path ordered products of 9 £ and %, for the Euclidean vacuum without 
the restriction that all the arguments are at the equal time, which will be a straight forward extension. In this manner, one can 
construct various excited states that are IR regular and free from the secular growth. 

B. Comparison to the recent publications 

In the recent papers ll38ll42ll . the absence of the secular growth is also claimed. It would be profitable to give a comparison 
between these works and our current work. First, in these papers, the item 1 raised in Sec. [I] i.e., the presence of the canonical 
system, which is equivalent to the original canonical system and whose interaction Hamiltonian is composed only of the IR 
irrelevant operators, is postulated, while this is not automatically guaranteed from the symmetry of the classical system. Second, 
in these papers, the mode function in de Sitter spacetime, whose amplitude is given by a constant Hubble parameter, is used in 
proving the conservation of the curvature perturbation. This leads to the quantitative discrepancy in the evaluation of the secular 

growth from ours. For instance, in Ref. ll38ll . the locality of the solution (x, t) given in Eq. (22) of the paper is crucial 
in their proof. However, the locality is not necessarily valid, once we take into account the fact that in the chaotic inflation, 
the amplitude of the fluctuation becomes larger and larger in the distant past as p oc e~ f dp£l . When we neglect this effect 
by setting (p/^/si) N in Eq. ( 14.23b to constant, the convergence condition is always satisfied (unless the interaction coefficient 
A, composed of the horizon flow functions, varies rapidly). Therefore our result does not contradict to the conservation of the 
curvature perturbation that they claimed. The third point is about the treatment of the UV contributions. In this paper, we 
have not directly discussed about the UV renormalization. We simply assumed that the UV divergent contributions, which are 
shown to be localized to the region where the two arguments of the Wightman functions are mutually almost light-like, can be 
renormalized by introducing the local counter terms. As long as the renormalization does not break the dilatation symmetry of 
the classical action, our discussion can hold. Recently, an interesting investigation about the UV renormalization is pursued in 
Ref. j42tl . It is claimed that a decaying composite operator in the free theory is kept decaying also after the renormalization 
of loops. Although the non-trivial assumptions such as the locality must be removed or verified, if this statement is correct, 
the conservation of the curvature perturbation can be shown also in the presence of the loop corrections. We should, however, 
emphasize that the conservation of the curvature perturbation does not prohibit the appearance of the logarithmic amplification, 
as we mentioned in Sec. |T] 

Finally, we also make a comment on the recent progress regarding the IR issues of a test field in the exact de Sitter spacetime, 
which can be interpreted as an approximation to the entropy mode. The regularity of the loop corrections for the Euclidean 
vacuum is shown for the massive scalar field by S. Hollands ll47ll and independently by D. Marolf and I. Morrison ll48l - l50ll . By 
contrast, for a massless scalar field, the IR regularity has not been shown and the absence of the secular growth is unclear Ipll - 
Isill (see also Ref. lf55tl ). Although the adiabatic curvature perturbation is a sort of massless field in the sense that the Wightman 
function G + (x, x') possesses the IR divergence and the long term correlation, the operation of the IR suppressing operators "R x , 
which appear by virtue of the residual gauge symmetry and by choosing the Euclidean vacuum, cures the singular behaviour. 
Hence, it would be intriguing to discuss a massless field with the exact shift symmetry in the de Sitter spacetime, in comparison 
with the case of the adiabatic mode. 
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Appendix A: Solving constraint equations 

In this section, we discuss about the boundary conditions of the constraint equations, which are elliptic type. By expanding 
the metric perturbations as £ = £j + £2 H — • > N = 1 + N\ + N2 H — • > an d = Ni l + iV» 2 H — • , the Hamiltonian constraint 
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and the momentum constraints yield 

VN n - 3/4 + e- 2 Pd 2 (n + pe- 2p d l N hn = H„ 
45 ?; (pN n - fn) ~ e~ 2p d 2 N hn + e^d^ 'Nj, r , 



Mi. 



(Al) 
(A2) 



where H n and Mi. n include n interaction picture fields (j in the combination (j — s or with differentiation. Eliminating N n 
from these constraint equations, we obtain 



1 - %- ) di&>N j>n - d 2 N^ n + 2^-e 2p d t ( n - 4^9i5 2 C = G i; „ , 



r 



v 



V 



where we defined 



4p. 



e 2p ( M i; „ - ^ff„ 



Operating 9* on Eq. ( IA3I ). we obtain 



We solve this equation as follows, 



d 2 8 l N hn = ^e 2p d 2 Cn - -d 4 ( n - ^d l C hn 



2p 2 



4p 2 



2p 2 



(A3) 



(A4) 



(A5) 



(A6) 



where G„(x) is an arbitrary solution of the Laplace equation, i.e., d 2 G^(x) — 0. Inserting this solution into Eq. (IA31) . we obtain 



d 2 N iin = d t 



A2 • 1 

n2 



2 . 2 e U ^ U ¥ 



2 (fl- 2 #c iin - G£) 



+ d t (n -;>*(■,„ - g£) - d. 



(A7) 



Again, introducing an arbitrary solution of the Laplace equation Gi. n (x), we solve Eq. dA7t as 



^e 2 "Cn(x) - -^d 2 Uz) ~ ^ 2 &C j<n {x) - G L n {x)) 



2p 2 y p ™ K ' 4p 2 

+ 9,3- 2 (d- 2 &C ]tn (x) - G L {x)) - d- 2 C hn (x) + G i>n (x) . 

Comparing the expression obtained by operating d l on Eq. dA8b with Eq. dA7t . we obtain 

fl'Gi.n = d l d- 2 C^ n - (d-z&C^n - G L ) . 
Using this expression, we rewrite the longitudinal part of G^„ as 

Gi, n = did' 2 [d l d- 2 C^ n - (d 2 iPC h „ - G L )\ + G itn - did- 2 d j G j<n . 
Inserting Eq. (lAlOt into Eq. (IA8b , we obtain 



(A8) 



(A9) 



(A10) 



N in (x) = did' 



2p 



2 e 2p Cn(x) - -a 2 C„(x) - (d~ 2 cPC j>n (x) - G L n {x)) 



4p 2 



(Si - did' 2 ®) {d- 2 C 3 , n {x) - Gj, n (x)) 



(All) 



When we perform quantization in the whole universe, it is natural to request the regularity of the perturbation at the spatial 
infinity. This requirement uniquely fixes G^ and the transverse part of Gi. n . Then, the shift vector depends on the curvature 
perturbation £ of the whole universe. To show the IR regularity, here we employ another boundary condition which requests 
that the integration region of the inverse Laplacian d~ 2 is confined to around the local observable region 0. As is shown in 
Refs. I25II26II . the degrees of freedom in changing the boundary condition can be understood as the gauge degrees of freedom 
in the local universe. Therefore, the operator 9 R is invariant under the change of the boundary condition. 
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Adjusting the solutions of the Laplace equations G^(x) and Gi. n (x), we can change the boundary condition for d 2 so that 



the integration region is limited. We fix the function G^ (x), requesting 

d- 2 W Lt {x)&C 3 , n {x) = d- 2 d j C j , n (x) - G%(x) . 



(A12) 



where we inserted the window function which takes a non-vanishing value only within the vicinity of the observable region 0. 
If we evaluate the term in the first line of Eq. (IA1 1) by using the Laplacian inverse with two different boundary conditions, 9-f 2 
and 2 , the difference is given by 



a 2 (d^ 2 [•••]-d l a 2 - 2 [•••]) =0, 



d l (flsSf 



did 2 



])=0. 



(A13) 



where we abbreviated the terms in the square bracket. Therefore, the change of the boundary condition for the Laplacian inverse 
can be absorbed by the transverse mode of Gi >n (x). Fixing the boundary condition of d~ 2 so that the integration region is 
restricted to the vicinity of the observable region, we obtain 



Ni, n {x) =d l d~ 2 W Lt {x) 



^ 2p Ux) - jd 2 Ux) - ^d- 2 W Lt (x)(PC j>n (x) 



(si - did- 2 ®) d- 2 w Lt (x)c^ n (x) . 



(A14) 



Inserting this solution into Eq. (IAU . we can also obtain the lapse function whose support of the Laplacian inverse d 2 is also 
confined. 



[1] D. Boyanovsky and H. J. de Vega, Phys. Rev. D 70, 063508 (2004) (arXiv: astro-p h/0406287 ] . 

[2] D. Boyanovsky, H. J. de Vega and N. G. Sanchez, Phys. Rev. D 71, 023509 (200 5) |arXiv:ast ro-ph/0409406 ]. 

[3] D. Boyanovsky, H. J. de Vega and N. G. Sanchez, Nucl. Phys. B 747, 25 (2006) |arXiv:astro-ph/0503669]. 

[4] D. Boyanovsky, H. J. de Vega and N. G. Sanchez, Phys. Rev. D72, 103006 (2005) |arXiv:astro-ph /0507596 ]. 

[5] N. C. Tsamis and R. P. Woodard, Annals Phys. 253, 1 (1997) |arXiv:hep-ph/9602316 1. 

[6] N. C. Tsamis and R. P. Woodard, Nucl. Phys. B 474, 235 (1996) |arXiv:hep-ph/9602315 |. 

[7] V. K. Onemli and R. P. Woodard, Class. Quant. Grav. 19, 4607 (2002) |arXiv"g r-qc/0204 d65]. 

[8] T. Brunier, V. K. Onemli and R. P. Woodard, Class. Quant. Grav. 22, 59 (2005) |arXiv:gr-qc/0408080 |. 

[9] T. Prokopec, N. C. Tsamis and R. P. Woodard, Annals Phys. 323, 132 4 (2008) larXiv:0707 0847i rgr-qcll. 
[10] M. S. Sloth, Nucl. Phys. B 748, 149 (2006) [arXiv:astro-ph/0604488|. 
[11] M. S. Sloth, Nucl. Phys. B 775, 78 (2007) |arXiv:h ep-th/06 12138], 
[12] D. Seery, JCAP 0711, 025 (2007) |arXiv:0707.3377 1 rastro-phlk 

[13] D. Seery, JCAP 0802, 006 (2008) larXiv:0707.3378 1 rastro-phll 

[14] Y. Urakawa and K. i. Maeda, Phys. Rev. D 78, 064004 (2008) |a rXiv:0801.0"l26l[hep-th]]. 

[15] P. Adshead, R. Easther and E. A. Lim, Phys. Rev. D 79, 063504 (2009) |arXiv:0809.4008 [hep-th]]. 

[16] H. R. S. Cogollo, Y. Rodriguez and C. A. Valenzuela-Toledo, JCAP 0808, 02 9 (20 08) I arXiv:080 6T546l rastro-phll . 

[17] Y. Rodriguez and C. A. Valenzuel a-Toledo, Phys. Re v. D 81, 023531 (2010) larXiv:08 1 1 .40921 rastro-phll . 

[18] D. Seery, JCAP 0905, 021 (2009) |arXiv:0903.2788l[astro-p h.CO]]. 

[19] X. Gao and F. Xu, JCAP 0907, 042 (2009) |arXlv :0905.0405 [hep-th] ]. 

[20] N. Bartolo, E. Dimastrogiovanni and A. Vallinotto, arXiv: 1006.0196 [astro-ph.CO]. 

[21] D. Seery, Class. Quant. Grav. 27, 124005 (2010) |arXiv: 1005. 1649 [astro-ph.CO]]. 

[22] E. O. Kahya, V. K. Onemli and R. P. Woodard, larXiv: 1 006. 39991 [astro-ph.CO] . 

[23] W. Xue, X. Gao and R. Brandenberger, JC AP 1206, 035 (2012) (arXiv: 1201.0768 [hep-th]]. 

[24] J. M. Maldacena, JHEP 0305, 013 (2003) |arXiv:astro- ph/02 1060 3], 

[25] Y. Urakawa, T. Tanaka, Phys. Rev. D82, 121301 (2010). larXiv: 1 007.0468 [hep-th]]. 

[26] Y. Urakawa and T. Tanaka, Prog. Theor. Phys . 125, 1067 (2011) larXiv: 1009.2947 1 rhep-thl 1 . 

[27] T. Tanaka and Y. Urakawa, larXiv:1209.1914i rhep-thl. 

[28] Y. Urakawa and T. Tanaka, Prog. Theor. Phys. 122, 779 (2009) |arXiv:0902.3209 [hep-th]]. 
[29] Y. Urakawa and T. Tanaka, Prog. Theor. Phys. 12 2, 1207 (2010) |arXiv:090 4.4415i rhep-thll. 

[30] Y. Urakawa, Prog. Theor. Phys. 126, 961 (2011) larXiv: 1105.10781 [he p-th] ] . 

[31] C. T. Byrnes, M. Gerstenlauer, A. Hebecker, S. Nurmi and G. Tasinato, larXiv: 1005.33071 rhep-thl . 
[32] M. Gerstenlauer, A. Hebecker a nd G. Tasinato,larXjv:1102.0560H astro-ph.COl. 
[33] S. B. Giddings and M. S. Sloth JarXiv: 1005. 105 6 1 [hep-th]. 

[34] S. B. Giddings and M. S. Sloth JarXiv: fT04 .0002 [hep-th]. 

[35] S. B. Giddings and M. S. Sloth, Phys. Rev. D 86, 83538 (2012) larXiv: 1 109.1000i rhep-thll. 
[36] L. Senatore and M. Zaldarriaga, larXiv: 1203.6354 [hep-th]. 

[37] G. L. Pimentel, L. Senatore and M. Zaldarriaga, JHEP 1207, 166 (2012) larXiv: 12 03.6651 [hep-th]]. 



23 



[38] L. Senatore and M. Zaldarriaga, arXiv: 12 10.6048 [hep-th]. 

[39] P. Creminelli, J. Norena and M. Simonovic, JCAP 1207, 052 (20 12) larXiv:1203.4595l [hep-th]]. 
[40] K. Hinterbichler, L. Hui and J. Khoury, JCAP 12 08, 017 (2012) larXiv:1203.6351i rhep-fhll. 
[41] T. Tanaka, Y. Urakawa, JCAP 1105, 01 4 (2011). larXiv:1103.123T1 rastro-ph.COll. 
[42] V. Assassi, D. Baumann and D. Green. larXiv: 1210.7792i rhep-thl. 

[43] A. Higuchi, D. Marolf and I. A. Morrison, P hys. Rev. D 83, 08 4029 (2011) larXiv:1012.3415i rgr-qcll. 
[44] D. Marolf, I. A. Morrison and M. Srednicki, larXiv: 1 209. 60391 rhep-thl . 
[45] N. C. Tsamis and R. P. Woodard, Annals Phys. 215, 96 (1992). 

[46] S. P. Miao a nd R. P. Woodard, JCAP 1207, 008 (2012) |arXiv:1204.1784i rastro-ph.COll. 

[47] S.Hollands, |arXiv:1010.53671 rgr-qcl. 

[48] D. Marolf and I. A. Morrison, Phys. Rev. D 82, 105032 (2010) larXiv: 1006.00351 [gr-qc]]. 
[49] D. Marolf and I. A. Morrison, Phys. Rev. D 84, 044040 (2011) |arXiv: 1010.5327 [gr-qc]]. 
[50] D. Marolf and I. A. Morrison, Gen. Rel. Grav. 43, 349 7 (201 1) |arXiv: 11044343] [gr-qc]]. 
[51] S. Hollands, Annales Henri Poincare 13, 1039 (2012) larXiv: 1105 . 1996 [gr-qc]]. 
[52] H. Kitamoto and Y. Kitazawa, Nucl. Phys. B 839, 552 (2010) larXiv: 1004.2451 [hep-th]]. 
[53] H. Kitamoto and Y. Kitazawa, Phys. Rev. D 83, 104043 (2011) |arXiv:1012.5930 [hep-th]]. 
[54] H. Kitamoto and Y. Kitazawa, Phys. Rev. D 85, 0440 62 (2012) larXiv: 1109 .4892 [hep-th]]. 
[55] A. Rajaraman, Phys. Rev. D 82, 123522 (2010) [arXiv: 1008. 1271 [hep-fh]]. 



